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Abstract — In this work we present two generators for the 
group of symmetries of the generic (m + 1) dimensional periodic 
Costas arrays over elementary abelian (Zp)™ groups: one that 
is defined by multiplication on m dimensions and the other by 
shear (addition) on m dimensions. Through exhaustive search 
we observe that these two generators characterize the group of 
symmetries for the examples we were able to compute. Following 
the results, we conjecture that these generators characterize the 
group of symmetries of the generic (m + 1) dimensional periodic 
Costas arrays over elementary abelian (Zp)™ groups. 



I. Introduction 

The original two dimensional Costas arrays were inttoduced 
in O, 131 and their periodicity properties have been studied 
in 0, in, and more recently in lITOl . Good correlation, and 
full periodicity properties such as those of some Costas arrays 
are very important for applications in digital watermarking lfT4ll . 
ifTSll . and optical communications Q, lfT3l . 

Multidimensional periodic arrays are especially useful for 
applications of digital watermarking for video, and combined 
video and audio media. In recent work Moreno et al. IfTTl 
presented a multi-periodical generalization of the Welch con- 
struction over the elementary abelian group (Zp)™. This 
Welch construction is too sparse for the application of digital 
watermarking, but in their work they show how its multi- 
periodicity property allows the substitution of columns by 
periodic sequences with good correlation properties. Three 
dimensional periodic arrays also have applications in optical 
communications. 

In this work, we generate many multidimensional periodic 
Costas arrays (MCPA) over elementary abelian (Zp)™ groups. 
These arrays are obtained through the application of alge- 
braic symmetries to the generalization of the multidimensional 
periodic Welch construction. We call these symmetries: Gi 
multiplication on the m dimensions symmetries and G2 shear 
(addition) on the m dimensions symmetries. We independently 
found Gi and G2 analysing the results obtained through ex- 
haustive exploration of some examples of these Welch arrays; 
later we found that these kinds of algebraic permutations were 
used on Quadratic residue arrays by Bomer et al [2], and on 
two dimensional Costas and sonar arrays [S]. 

We compared the number of results obtained from the 
exhaustive exploration of the examples (Za)^, (25)^, and (22)"^ 
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with the number of results obtained using these generators; 
and observed that the group of symmetries for the examples is 
characterized with the use of the algebraic generators. In other 
words, all the MPCA found through exhaustive exploration are 
generated using a combination of the symmetries that will be 
presented in this work. 

The rest of the paper is organized as follows: Section II- A I 
contains the formal definitions of a generic (m+l)-dimensional 
periodic Costas array over the elementary abelian group (Zp)™, 
the definition of the Costas property, and reviews two algebraic 
symmetries of multidimensional periodic Welch Costas arrays 
that are used in combination with the symmetries presented 
in this work to characterize the group of symmetries of the 
generic MPCA. Sections HIl and HVl present the generators used 
to complete the characterization of the group of symmetries of 
the generic 3D MPCA. Sections |lll| and [V] present the gener- 
alization of the generators for multiple dimensions. Section IVll 
describes the computational work performed; and section IVIII 
concludes this work, with our conjectures on the existence of 
the generic multidimensional periodic Costas arrays. 

A. Definitions 

We begin by providing the definitions of the generic and 
general (to + l)-dimensional periodic Costas array over the 
elementary abelian group and introduce the Generalized Welch 
construction. 

Definition 1. A generic (to + 1) -dimensional periodic Costas 
array over the elementary abelian group (Zp)™ is a permutation 
function f : ((Zp)™)* (Zpm_x), where A* means A-{0}. 
This function has the distinct difference property.- for any h ^ 

0, a, 6 e (Zp)", f{a+h)-f{a) ^ f{b+h)^f{b) implies a ^ b, 
where the addition and subtraction operations are performed 
in the corresponding abelian group. 

Example 1. The following is a grid defined over Z5 x Zs.- 
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The following is a three dimensional periodic Welch Costas 
array over the elementary abelian group Z5 x Z5.- 
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Remark: Through the examples of this work the rows 
are numbered upwards, starting from the bottom, whilst the 
columns are numbered from left to right. 

Definition 2. A general (m+1) periodic Costas array is defined 
as before, but * is now over any grid location. 

The well known Welch construction is generalized to con- 
struct MPCA. 

Construction 1. (Generalized Welch: a multi-periodic Costas 
array in m+1 dimensions with correlation 1 ). 

Consider the sequence s such that s^i — loga{X) where 
X G GF(jj"^), and a is a primitive element of GF(jj"^). 
Specifically for X = a*, Sq,; = logaoi^ — i. The sequence 
scheme is to use the m-tuple representation of a' to determine 
the coordinates (location) on the m dimensional integer grid 
defined above, s is a periodic sequence with period — 1. 
For our Costas type construction we take the grid point location 
belonging to s^i, and place a 1 in a column of length 1 
located above the grid point, and zeros in all other entries in 
that column. 

Tlieorem 1. The permutation s posseses the Costas property: 
S(a+h) — Sa — S(6+h) — Sfc — fl = 6, V/i ^ and is therefore 
a Costas array. 

Proof: 

Let a,b,he GFip'"^) and /i 7^ 



Theorem 3. Addition ( modulo — 1 j 0/ any integer less than 
p"^ — \ to a periodic Costas array generates a new periodic 
Costas array. 

Example 3. Addition of 4 to W. 
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The addition and multiplication (modulo — 1) symmetries 
are used in combination with the symmetries presented in the 
next sections to generate many MPCA and to characterize the 
group of symmetries of the examples that we also characterized 
computationally. 

II. Row AND Column Multiplication 

For the 3D Welch MPCA, the rows and columns multiplica- 
tion symmetries consist in shifting the rows and columns of the 
two dimensional representation. In general, these symmetries 
consists of shifting m-dimensions of the (m + 1) dimension 
Welch MPCA x times their original position. In the following 
sections we define the row multiplication (RGl), column mul- 
tiplication (CGI) and the combination of the row and column 
multiplication symmetries (Gl), and apply them to the 3D 
Welch MPCA. Finally we define a generalization of the Gl 
symmetry for (m + 1) dimensions. 



log{a + h) — log{a) = log(b + h) — log(b) (1) 

h h 

! + - = ! + - (3) 

a 

Therefore, for h ^ Q ^ a = b. 



B. Addition and Multiplication ( modulo p™ — 1) symmetries 

The following two algebraic symmetries were introduced by 
Moreno et al. in a previous work [12] and are used to character- 
ize the group of symmetries of our examples. The multiplication 
(modulo — 1) is the first symmetry applied to the sequences 
obtained from the Welch construction, then we apply the Gl 
and G2 symmetries introduced in this work, and finally we 
use the addition symmetry to finish the characterization of our 
examples. 

Theorem 2. Multiplication (modulo p™ — 1) of a periodic 
Costas array by an integer less than and relatively prime to 
p™ — 1 generates a new periodic Costas array. 

Example 2. Multiplication of W by 23 = -1 mod 24. 

/ 23 7 10 9 14 \ 

17 3 1 8 4 

5 16 20 13 15 

11 2 21 22 19 

\ * 18 6 12 / 



A. Row Multiplication 

Let G (Zp)^, and x E {1 . . .p — 1}, and let (w^.j) be 
the element of matrix W in position 

Definition 3. (Row Multiplication) 

RGlx '. (^p)^ ^ (^p)^ ^ permutation function such that 
RGlx '■ (wi.j) — ^ {Wxi,j) permutes the rows from position i to 
position xi of a three dimensional periodic Costas W over the 
elementary abelian group Zp x Zp. 

Example 4. The following three dimensional periodic Welch 
Costas array over the elementary abelian group Z5 x Z5 is 
one symmetry, out of 4 distinct row multiplications symmetries 
obtained after applying RGl2{W). 



RGhiW) = 
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B. Column Multiplication 

Let G (Zp)^, and x G {I . . .p — 1}, and let (w^j) be 
the element of matrix W in position 

Definition 4. (Column Multiplication) 

CGlx '. (^p)^ — ^ (^p)^ "2 permutation function such that 
CGlx ■ (wij) — > {wi_xj) permutes the columns from position j 
to position xj of a three dimensional periodic Costas W over 
the elementary abelian 
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Example 5. The following three dimensional periodic Welch 
Castas array over the elementary abelian group Z5 x Z5 is one 
symmetry, out of 4 distinct column multiplications symmetries 
obtained after applying the CGl3{W). 



RGhiW) = 
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We can combine both transformations to generate more 
Welch Costas arrays. 

Note that RG1,,(CG1,,(W)) = CGI,, (RGl,, (W)). 

Example 6. The following is the three dimensional periodic 
Welch Costas array over the elementary abelian group Z5 x Z5 
after applying the RGl2{CGlz{W)). 



RGUiCGhiW)) 



C. Row and Column Multiplication 

Definition 5. (Row and Column Multiplication) Glxi,x2 
(Zp)^ — 5> (Zp)^ is a permutation function such that Glxi,x2 ■ 
~^ {^xii,x2j) permutes the rows from position i to 
position xii and the columns from position j to position X2j 
of a three dimensional periodic Costas W over the elementary 
abelian group Zp x Zp. 

Example 7. The following is the three dimensional periodic 
Welch Costas array over the elementary abelian group Z5 x Z5 
after applying the Gl2,3(W). 
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Theorem 4. The application of the permutation function 
Gl xi,x2 to a three dimensional periodic Costas array over 
the elementary abelian group Zp x Zp generates a new three 
dimensional periodic Costas array over the elementary abelian 
group lip X Zp. 

Proof Let {i' {hi, /12) G Zp x Zp, and Xi,X2 G 

{0...P-1}, (a;i,a;2)/(0,0) 



log(xiia + X2j + hia + /i2) — log{xiia + X2j) 



log 



xiia + X2j + hia + 
x\ia + X22 



log\\^ 



hia + h2 
xiia + X2j 
Proof of the Costas property 



(4) 



(5) 



(6) 



logll 



hia + /i2 



xiia 
log (^1 -f 

hia + 



■X2j 

hia 



xii'a + X2j' 
h2 



(7) 



xiia + X23 
hia 



1 + 

xii'a + X23' 
Therefore for {h^ h2) + (0,0) (1,3) = 



III. Generalization: (to + 1)-dimensions 

MULTIPLICATION PERMUTATION 



(8) 



Let a 



(oi, 02, . . . , a„i) 



and 



{xi,X2, ■ . ■ ,Xm) G (Zp)'", such that Xi e {l...p— 1}, 
and let (Wai,a2,....a„) be the elements of m dimensional 
matrix W in positions (ai, 02, ... , a„i). 



Definition 6. (Gl: Multiplication Permutation) 

Glxi,x2,...,Xm-{'^pY^ ^ (^p)™ '■^ permutation function such 

that Gl xi ^X2,. ■ ■ ,Xjn- i}^ ai ,a2 . .,ajn^ ^ (}^xiai,X2a2....,Xjnajn^ per- 
mutes the elements of matrix W from position (oi, 02, ... , flm) 
to position (xiai, X2CL2, ■ ■ ■ , XmO-m) of a (to + 1) dimensional 
periodic Costas W over the elementary abelian group (Zp)"'. 

Theorem 5. The application of the permutation function Gl 
to a (m + 1) dimensional periodic Costas array over the 
elementary abelian group (Zp)™ generates new (m + 1) di- 
mensional periodic Costas arrays over the elementary abelian 
group (Zp)™ 



Proof: Let h,a,b € GF{p"') and x e (Zp) 
G {1 . . .p—1} and let (x-a) — (xiuia™- 



-X2a2a 



such that 

m-2 I 



log 



log{{a ■ x) + h) — log{x ■ a) 
(x ■ a) + h\ , h 



X ■ a 

Proof of the Costas property 



= log[l 
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log [I 



[x ■ a) 
h 



= log 1 
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Therefore for h ^i) ^ a = h 



= 1 
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IV. Row AND Column Shear 

The row (column) shear symmetry consists shifting each 
of the elements of the row (column) of the two dimensional 
representation of the 3D Welch MPCA by a different number 
of shifts. Example |8] shows how to obtain a new MPCA by 
shifting rows through 4 by through 4 times to the right (E.g 
row 1 is cyclically shifted 1 time to the right, row 2 is cyclically 
shifted 2 times to the right, and so on). In the following sections 
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we define the row shear (RG2), column shear (CG2), and the 
combination of the row and column shear symmetries (G2), 
and apply them to the 3D Welch MPCA. Then we define a 
generalization of the G2 symmetry for (m +1) dimensions. 

A. Row Shear 

Let (i, j) G (Zp)^, and x G {0 . . .p — 1}, and let (wjj) be 
the element of matrix W in position 

Definition 7. (Row Shear) 

RG2x '■ (Zp)^ — > (^p)^ is a permutation function such that 
■ ~^ ^i,j+xi shifts the elements of the rows i from 
position j to position j + xi of a three dimensional periodic 
Costas W over the elementary abelian group "Lp x Zp. 

Example 8. The following three dimensional periodic Welch 
Costas array over the elementary abelian group Z5 x Z5, is 
one symmetry, out of 5 distinct column stairlike shift symmetries 
obtained after applying the RG2i(W). 



RG2i{W) 



B. Column Shear 

Let (i, j) e {"LpY, and .x G {0 . . .p - 1}, and let (Wjj) be 
the element of matrix W in position 

Definition 8. f Column Shear) 

CG2x : (Zp)^ (^p)^ be a permutation function such that 
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CG2, 



w 



i+xj,j 



shifts the elements of the columns 



j from position i to position i + xj of a three dimensional 
periodic Costas W over the elementary abelian group 

Example 9. The following three dimensional periodic Welch 

Costas array over the elementary abelian group Z5 x Z5, is 
one symmetry, out of 5 distinct column stairlike shift symmetries 
after applying the CG22{W). 



CG22{W) = 



Note that because in the permutation function, the new 
position in the first dimension depends on the position of 
the second dimension RG2^(CG2y(W)) ^ CG2j,(RG2^(W)). 
However the composition of these generators generates another 
MPCA. 

Example 10. The following is the three dimensional periodic 
Welch Costas array over the elementary abelian group Z5 x Z5 
after applying RG2i {CG22{W)). 
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Let G Zp X Zp, and x\,X2 G {0. ..p — 1}, and let 
(Wjj) be the element of matrix W in position 

Definition 9. (Row and Column Shear) G2xi^x2 '■ (^p)^ 
(Zp)^ is a permutation function such that G2x-^.X2 ■ i^i.j) 
(wi-^x2j,j+xii) shifts the elements from position to posi- 
tion (i + X2j,j + xii) of a three dimensional periodic Costas 
W over the elementary abelian group Zp x Zp whenever 
{'i + X2j,j + xii) 7^ {0}. 

Theorem 6. The application of the row and column shift 
permutation function G2xi,x2 to a three dimensional periodic 

Costas array over the elementary abelian group "Lp x Zp 
generates a new three dimensional periodic Costas array over 
the elementary abelian group 

Proof: Let (i, j), {hi, /i2) G Zp x Zp, and xi,X2 G 

{0...P-1} 



log{{i + X2j)a + j + xii + hia + /12) 

—log{{i + X2j)a + J + xii) = 
{i + X2j)a + j + xii + hia + /12 



log 



log [I 
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hia + /?,2 
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Proof of the Costas property 
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hia + ft,2 



(i' + X2j')a + j' + xii' 



{i + X2j)a + j + xii 

= {i' + X2f)a + j' + xii' 
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a{i - i') + ax2{j - f) +xiii- i') + (j - /) = (19) 
Therefore for {hi,h2) + (0,0) ^ (z, j) = (z', j') 



V. Generalization: (m + 1)-dimensions shear 

PERMUTATION 

Let a = (fli, 02, . . . , am) G (Zp)™, fc, r G {1 . . . m}, k ^ r, 
and X G {1 . . — 1}, and let (Wai,a2,...,a,„) be the elements 
of m dimensional matrix W in positions (ai, 02, . . . , am). 

Definition 10. {G2: Shear Permutation) 

G2:{Zp)"^ (^p)™ is d permutation function such that 
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permutes the elements of matrix W from position 
(ai, 02, . . . , Cm) to position (ai + 02^2 + •■• flmUm, OlWl + 
12 + • • ■ CLmVm + • . . ) of a {ni + 1) dimensional periodic 
Costas W over the elementary abelian group (Zp)"' whenever 

(ai + a2M2 + . . . flmUm, fll^l + 02 + . . . a„iVra + . . . ) 7^ {0}- 

Theorem 7. The application of the permutation function G2 
to a (m + 1) dimensional periodic Costas array over the 
elementary abelian group (Zp)™ generates new (m + 1) di- 
mensional periodic Costas arrays over the elementary abelian 
group (Zp)™ 

Proof: The proof is similar to the proof in Theorem |6] ■ 

VI. Computational results 

The exhaustive exploration for the enumeration of multidi- 
mensional periodic Costas arrays (MPCA) has the same compu- 
tational complexity of the well known problem of enumeration 
of two-dimensional Costas arrays of size — 1, which is 
(p™ — I)!- More precisely the exhaustive exploration for the 
enumeration of periodic MPCA over the elementary abelian 
group is (p™ — 2)!; because we can fix the first defined (i.e. 
non zero position) to a constant value k and then permute all the 
other positions. This is due to the fact that the arrays are multi- 
periodic and therefore all multi-periodic shifts are equivalent. 
All the other arrays with values different tofc are obtained by 
the addition symmetry. 

At the moment of writing this paper, the last enumeration 
reported finished for two dimensional Costas arrays is for 
length 28 14| and the time per single CPU was determined 
to be 70 years. With our resources we were able to enumerate 
the multidimensional periodic Costas arrays over elementary 
abelian groups of length (Zs)^, (Zs)^ and (Z2)^, which are 
similar to computing the two dimensional Costas arrays of 
lengths 7, 23 and 6 respectively. The complete enumeration 
for (Zy)^ or (Za)^ would require approximately the equivalent 
to compute the enumeration of two dimensional Costas arrays 
of length 47 and 25 which will require high performance 
computing to be completed in a reasonable amount of time. 

A. The algorithm 

The algorithm to compute the multidimensional examples is 
similar to the algorithm to construct two dimensional Costas 
arrays |[T], (|4]. It consists in a one dimensional permutator 
that uses a backtracking algorithm to prune the search space. 
Backtracking is a general algorithm for solving computational 
problems by incrementally generating all possible solutions. 
The execution of a backtracking algorithm can be modelled 
as a search tree where every node is a partial solution. Moving 
forward corresponds to approaching a valid solution, and going 
backwards corresponds to abandoning a partial candidate that 
cannot possibly generate a valid solution. In this case the 
possible solutions are the MPCAs, and the partial solutions are 
sub permutations that meet the multi-periodic Costas property. 
The main change to the original algorithm for two dimensions 
is in the function to check whether a sub permutation meets 
the three dimensional or four dimensional Costas property. We 
have parallel implementations of these algorithms using MPI, 
CUDA (CPUs), and FPGAs. 



VII. Conclusions 

We have presented two generators Gl and G2 for the 
group of symmetries of MPCAs over elementary abeUan group 
(Zp)™. By applying Gl and G2 we characterized the full 
group of symmetries for the examples (Za)^, (Zs)^ and (Z2)'^. 
Because of the strong properties of MPCA, perfect autocor- 
relation and multi-periodicity, we conjecture that generators 
Gl and G2 characterize the full group of symmetries of the 
multidimensional periodic Welch Costas arrays over elementary 
abelian group (Zp)™. 

We also completed exhaustive search for the cases: (Z4)^, 
(Z2 X Z3), (Z3 X Z4) with empty results, and for the same reason 
as before, we conjecture that there are no multidimensional 
periodic Costas arrays other than multidimensional periodic 
Costas arrays over the elementary abelian group (Zp)™. 
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